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Abstract. We present here the results of a detailed study of the X-ray properties of the cluster of galaxies Abell 3128 (z = 0.06), 
based on the analysis of deep (100 ks) XMM-Newton data. The most obvious feature of the X-ray morphology of A3128 is the 
presence of two X-ray peaks separated by ~12'. By detecting the redshifted Fe K line, we find that the Northeast (NE) X-ray 
peak observed toward A3 128 is a distant luminous cluster of galaxies at redshift z = 0.44. Our subsequent optical spectroscopic 
observation of a distant radio bright galaxy in the centre of the NE X-ray peak with the Magellan telescope also revealed a 
redshift of z = 0.44, confirming the association of the galaxy with the cluster seen in X-rays. We detect a gravitational arc 
around the galaxy. The properties of this galaxy indicate that it is the cD galaxy of the cluster in the background. The properties 
of the Southwest X-ray peak suggest that it is the core of a group merging with A3 128 along our line of sight. Based on 2D 
maps of thermodynamic properties of the intra-cluster medium determined after subtracting a model for the background cluster, 
we conclude that an enhanced surface brightness region at a distance of ~2.8' from the centre of the galaxy distribution is the 
centre of the gravitational potential of the cluster A3128. The unrelaxed nature of A3128 can be attributed to its location in the 
high density environment of the Horologium-Reticulum supercluster. 
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1. Introduction 

The formation and growth of the largest structures in the 
Universe, clusters of galaxies and superclusters, is a still ongo- 
ing process. According to the standard cosmological scenario 
clusters of galaxies continuously form and grow through merg- 
ing with groups and individual galaxies, which are falling in 
along filaments. In several cases clusters are found close to 
each other in the highest overdensity regions of the Universe: 
the superclusters. In the crowded environment of superclusters, 
we can often study processes related to the growth of structure, 
like mutual interaction between the member clusters and their 
accretion history. 

The Horologium-Reticulum (H-R) supercluster (z — 0.06), 
together with the better known Shapley Supercluster are the 
most massive structures in the local universe within a dis- 
tance of 300 Mpc, with an estimated total mass of 10 17 M . 
Einast o et al. (2001) list 35 member clusters in H-R. 

Abell 3 128 is a rich, highly substructured cluster in the H-R 
supercluster. It has been well studied using redshift and radio 



surveys, and by Ch andra (20 ks: |Rose et al. 1 |2002|) . Based on 
A^-body simulations. ICaldwell & Rosd d 19971) proposed that in 
the past A3 128 encountered a tidal interaction with the clus- 
ter A3125 (current separation 1°, corresponding to ~6 Mpc). 
Rosat and Chandra observations of A3 128 show a complex, 
disturbed X-ray morphology. Embedded in a ~20' diameter dif- 
fuse halo there are two cores in A3 128, separated by 12'. The 
Northeast (NE) and Southwest (SW) cores comprise ~85% of 
the X-ray emission of A3 128, the remainder being in a diffuse 
halo. While the more luminous SW core is centred on a com- 
pact group around a bright elliptical galaxy, the NE core does 
not coincide with any bright galaxy. The NE core is asymmetric 
and elongated, with the peak of the emission slightly displaced 
to the Northeast. With core radii of ~30 kpc, both cores are un- 
usually narrow as compared to other non-cooling core clusters 
with similar temperatures (~3.5 keV). 

Ros e et al. I (l2002l) proposed a model to explain the double- 
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peaked nature of the X-ray emission. In position-position and 
position-redshift diagrams they identified an infalling group 
and a candidate for a post-passage tidally distended group (fil- 
ament), which they propose as candidates for producing the 
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complex structure seen in X-ray images of the intra-cluster 
medium (ICM). They explain the SW core as the still intact hot 
ga s of an infalling group. The NE core represents according 
to iRose et al.1 ( 12002b the surviving ICM of a group of galax- 



ies that has fallen supersonically (Mach number 6) into the 
cluster A3 128 along the main filament connecting A3 128 and 
A3 125. The infalling group has, according to this scenario, 
passed through the core of A3 128 and the galaxies of the group, 
moving ballistically through the cluster, are by now well ahead 
of the stripped gas. The high infall velocity is caused by the 
deep potential well of the H-R supercluster, which makes even 
an infall of a small group an energetic and inter esting event 
in the life of a galaxy cluster. IRose et al.1 (12002b also report 
a possible detection of a radio arc, which lies slightly to the 
Northeast of the peak emission in the NE core, where the bow 
shock is expected if the gas is moving supersonically. 

Howev e r, the Fe abundance of the NE core found by 
Rose et al.1 (120021) is low (0.13 solar) and barring conclusive 
redshift information its association with A3 128 could not be 
confirmed. As an alternative scenario to their u nified view 
of the merging events occurring in A3128/A3125. lRose et al 



(2002) invoke the possibility that the NE emission peak is as- 
sociated with a background cluster and the detected radio emis- 
sion is associated with a radio bright central galaxy in a back- 
ground cluster. In order to obtain a conclusive measure for plac- 
ing of the NE component, a deeper observation for a redshift 
measurement using X-ray emission lines was needed. 

Here we present the results of a deep 100 ks observation 
of Abell 3128 with the European Photon Imaging Cameras 
(EPIC) on XMM-Newton. The observation allows us to mea- 
sure the redshift of the NE core and to unambiguously demon- 
strate that the dominant fraction of its emission originates in a 
background cluster. We also present here data obtained by the 
Magellan telescope, that allow us to confirm the association of 
the X-ray emission with a background cluster surrounding a ra- 
dio bright cD galaxy. The large effective area of XMM-Newton 
combined with the deep exposure, provides us with sufficient 
statistics to map the 2D distribution of thermodynamic proper- 
ties of the ICM and thus to study the merging history of the dis- 
turbed, complex cluster of galaxies A3 128. The XMM-Newton 
Reflection Grating Spectrometer data have insufficient statis- 
tics to derive accurate spectral properties for the cluster. 

Throughout the paper we use Hq = 70 kms^'Mpc , 
Qm = 0.3, £2a = 0.7, which imply a line ar scale of 78 kp c 



at the redshift of A3 128 (z = 0.06, IRose et al.ll2 002) 



Unless specified otherwise, all errors are at the 68% confidence 
level for one interesting parameter (A^ 2 = 1). The elemental 
abundances are given with respect to the proto-solar values of 
Loddersl(l2003l) . 



2. Observations and data analysis 

2.1. XMM-Newton data 

Abell 3128 was observed with XMM-Newton during two 
pointings on May 29-30 and between May 31 and June 1st 
2006 (revolutions 1185 and 1186) for a total exposure time of 
~100 ks. The two pointings were centred on the NE (72 ks) 



and on the SW core (32 ks), respectively. The EPIC/MOS de- 
tectors were operated in the full frame mode, while for the 
EPIC/pn detector the extended full frame window mode was 
employed. The observations were performed using the thin fil- 
ter. The calibrated event files were obtained using the 7.0.0 ver- 
sion of the XMM-Newton Science Analysis System (SAS). For 
EPIC/MOS we keep only the single, double, triple, and quadru- 
ple pixel events (PATTERN<12), while for EPIC/pn, we make 
use of singe and double events (PATTERN <4). The spectral 
redistribution and ancillary response files are created with the 
SAS tasks rmfgen and arfgen separately for each camera and 
spectral extraction region that we analyze. 

Because of low number of counts in original bins, the ex- 
tracted spectra are rebinned into bins with a minimum of 30 
counts per bin. Our rebinned spectrum has bins of size at least 
1/4 times the FWHM of the instrument. We fit the MOS1, 
MOS2, and pn spectra from both pointings simultaneously with 
the same model, with their relative normalizations left as free 
parameters. Our spectral analysis is restricted to the 0.4-7.0 
keV band. In the spectral analysis we remove all bright point 
sources with a flux higher than 4.8 x 10~ 14 ergs s" 1 cm" 2 . 

For the spectra l analysis we use the SPEX package 
(IKaastra et al.l 1 1996b . We fix the Galactic absorption in our 
model to t he value deduced from the Hi data (Nu = 1.47 x 
10 cm , Dickey & Lockmaniri 990). We fit the cluster spec- 
tra with a plasma model in collisionally ionized equilibrium 
(MEKAL). The free parameters in the MEKAL model are the 
normalisation, temperature, and the Fe abundance. Since for 
elements other than Fe we can not obtain accurate abundance 
values, we fix them to 0.4 solar in the model. 



2.2. Background modeling 

2.2.1. Particle and instrumental background 

To clean the data from the soft proton induced events, we ex- 
tract light curves for each detector separately in the 10-12 keV 
energy band where the cluster emission is negligible and the 
detected emission is dominated by the particle induced events. 
Since flares having a particularly soft spectrum may be missed 
when only the high energy part is studied, we also extract 
lightcurves in the 0.3-2.0 keV band. A visual inspection of the 
light curves reveals that the observation is not badly affected by 
the soft proton flares. After excluding the time periods when the 
count rate in the two considered energy ranges deviates from 
the mean by more than 3<x, we are left for the first pointing 
with 7 1 ks and 67 ks, and for the second pointing with 3 1 ks 
and 25 ks for EPIC/MOS and EPIC/pn, respectively. 

We subtract the EPIC instrumental background using 
closed-filter observations. For EPIC/MOS we use a closed fil- 
ter observation (obs. ID 0150390101) with an exposure time of 
200 ks. For EPIC/pn we use a closed filter observation (obs. ID 
0106660401) with an exposure time of 120 ks. The instrumen- 
tal background consists of fluorescence line emission, intrinsic 
instrumental noise, and particle induced noise caused by high- 
energy cosmic rays which are able to reach the detector even 
when the filter wheel is in closed position. The instrumental 
background varies from observation to observation. In order to 
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Table 1. The CXB components used in the fitting of the A3 128 
spectra. The power-law photon index is F = 1.41. The unab- 
sorbed fluxes are determined in the 0.3-10.0 keV band. 



Comp. 


kr 


Flux 




(keV) 


(10" 12 erg cm" 2 s" 1 deg" 2 ) 


LHB/SDC 


0.08 


2.5 ± 0.6 


HDC 


0.25 ± 0.02 


4.0 ± 0.6 


EPL 




22.0 ± 1.9 



scale the closed filter observations to the instrumental back- 
ground during the source observation, we use the events regis- 
tered outside the field of view (out-of-FOV) of the EPIC detec- 
tors, outside a radius of 15.4' from the FOV centre. Separately 
for each instrument, we scale the closed filter observation by 
the ratio of the 8-12 keV out-of-FOV count rate in our cluster 
observation and in the closed filter observation. 

2.2.2. The cosmic X-ray background 

We correct for the Cosmic X-ray Backgro und (CXB) during 
spectral fitting. Kuntz & Snowden (2000) distinguish 4 dif- 
ferent background/foreground components: the extragalactic 
power-law (EPL), the local hot bubble (LHB), the soft dis- 
tant component (SDC), and the hard distant component (HDC). 
The EPL component is the integrated emission of faint dis- 
crete sources, mainly distant Active Galactic Nuclei (AGNs). 
The LHB is a local supernova remnant, in which our Solar 
System resides. It produces virtually unabsorbed emission at 
a temperature of ~10 6 K. The soft and hard distant components 
originate at larger distances. They might be identified with the 
Galactic halo, Galactic corona or the Local group emission and 
are absorbed by almost the full Galactic column density. Using 
the spectral band above 0.4 keV we can not reliably distin- 
guish the SDC emission from the LHB component. Therefore, 
at temperatures below 0. 1 keV we only consider the contribu- 
tion of one thermal component. 

The large field of view, with regions where the contribu- 
tion of cluster emission is small, enables us to determine the 
local properties of the background emission. We fit the spec- 
tra extracted from a region where the contribution of cluster 
emission is small and the emission is dominated by the X- 
ray background. We model the EPL emission with a power- 
law with a photon index of 1.41 dDe Luca & Molendil 120041) . 
We model the soft foreground emission (LHB/SDC and HDC) 
by 2 collisionally ionized plasma models (MEKAL). Since 
the fitted spectral band does not allow us to accurately con- 
strain the LHB/SDC temperature, w e fix its value to 0.08 keV 
(based on iKuntz & Snowdenll2000l) . We leave the HDC tem- 
perature as a free parameter in the fit. To account for the re- 
maining cluster emission in the extraction region, the contri- 
bution of which to the total flux is ~30%, we use an addi- 
tional thermal model. We find that both the temperature and 
the flux of the emission attributed to the HDC (0.25 keV 



data (0.127 keV and 2.44 x 10" 12 erg s" 1 cm" 2 deg" 2 , respec- 
tively). This difference might be due to the contribution of the 
intra-supercluster medium. The adopted values of the tempera- 
tures and of the 0.3-10.0 keV fluxes of the background compo- 
nents are given in TableQ] In our subsequent spectral fits we fix 
the parameters of the background components to these adopted 
values. 



2.3. The X-ray images 

The spatial resolution of XMM-Newton allows a robust image 
reconstruction down to 8" scales. However, due to the large 
wings of the point spread function (PSF), there is a significant 
contamination from the emission on scales of 8" or smaller 
to large spatial scales upto several arcminutes. For the pur- 
pose of cluster survey work (Finoguenov et al. in prep.), we 
have developed a procedure, which performs an image restora- 
tion using a symmetr i c mod el for the XMM PSF and the cal- 
ibration of iGhizzar di (2001). The flux on small scales is esti- 
m ated performing t he sc ale-wise wavelet analysis, as described 



and 4.0 x 10" 12 erg s" 1 



cm 2 deg z , respectively) are higher 



than the valu es determine d based on the results published by 
IKuntz & Snowdenl d2000h using the ROSAT All Sky Survey 



in IVikhlinin et al.l (119981) . Then the estimated flux is used to 
subtract the PSF model prediction on large scales and increase 
accordingly the flux on the small scales. The small scales in 
this procedure are a sum of 8" and 16" scales, which allows 
to avoid the variation in the PSF shap e with off-axis angle, as 
described in iFinoguenov et al.l (120071 COSMOS special issue 
paper). 

In Fig.l we display both the large scales of the emission 
starting at the 32" scale together with smaller scales (8" and 
16"). 

3. Results 

The most obvious feature of the X-ray morphology of A3 128 
is the presence of two X-ray peaks separated by ~12'. In Fig.Q] 
we show the X-ray contours superimposed on the Digitalized 
Sky Survey image of the cluster. The thick contours highlight 
the large scale (32") features in X-ray morphology, the thin 
contours highlight the features on small spatial scales (8" and 
16"). 

While the SW peak of the X-ray intensity coincides with 
a bright elliptical galaxy (ENACS 75) surrounded by a com- 
pact group of galaxies, the NE core does not coincide with any 
bright galaxy. On the small spatial scales the NE peak has a 
strongly elongated morphology and on larger scales it seems 
to have a tail toward the SW direction. The thick contours also 
reveal an association of the hot gas with two small groups of 
galaxies to the East of the main cluster contours. 

3.1. Properties of the X-ray peaks 
3.1 .1 . The global properties 

In order to study the global spectral properties of the two X-ray 
bright peaks, we extract their spectra from circular regions with 
a radius of 1.5'. In Fig. |2]we show the spectra with their best 
fit model. Previous Chandra data showed that while the SW X- 
ray core has a high Fe abundance, the Fe abundance of the NE 
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Fig. 1. X-ray contours superimposed on the Digitalized Sky Survey image of A3 128. North is up, West is to the right. The thick 
contours highlight the large scale features in the X-ray morphology, the thin contours highlight the features at small spatial scales. 
The contour spacing is arbitrary. 



core is very low, and therefore it is also not possible to confirm 
its association with the cluster A3 128 using the redshift of the 
Fe K line. Even using the deep XMM-Newton data we do not 
see a clear Fe K emission line at the expected energy for the 
redshift of A3 128. However, we detect in all EPIC detectors 
for both XMM-Newton pointings a line at ~4.6 keV. This cor- 
responds to the Fe K line emission of a cluster at a redshift of 
-0.45. 

We also extract a spectrum from a circular region with a 
radius of 2' centred on the apparent tail of the NE X-ray peak 
directed to the SW (a = 3 h 30 m 30 s , 5 = -52°31'0"). We im- 
mediately see an Fe K line at ~6.3 keV and we see no obvious 
line feature around 4.6 keV. We fit the spectrum extracted from 
the "tail" region and from the SW X-ray peak with a thermal 
model. The free parameters of the fit are the redshift, tempera- 
ture, Fe abundance, and emission measure. The best fit param- 
eters are shown in Table [2] 

We fit the spectrum of the NE peak with a combination of 
two thermal plasma models: one for the emission of A3 128 
and the other for the emission of the background cluster. We 
fix the parameters of the model of the foreground cluster to the 



values determined for the tail region: redshift z — 0.06, metal- 
licity Z = 0.49, and temperature kT = 3.5 keV. The emission 
measure of the foreground cluster and the redshift, Fe abun- 
dance, temperature, and emission measure of the background 
cluster are free parameters in the spectral fit. The best fit red- 
shift value for the background cluster is z = 0.444 ± 0.003. 
The linear scale at this redshift for the ACDM cosmology is 
341 kpc arcmin -1 . The best fit temperature, Fe abundance, and 
emission measure of the background cluster are shown in the 
column "NE" of Table [2] Our spectral analysis indicates that 
~60% of the X-ray emission in the extraction region is coming 
from the background cluster. The 0.3-10.0 keV luminosity of 
the background cluster within the radius of 1.5' (512 kpc) is 
4.7 x lO^ergss" 1 . 

The best fit parameters of the background cluster depend on 
the assumptions about the emission of the foreground cluster. 
Therefore, the systematic uncertainties on the best fit param- 
eters are larger than the quoted statistical uncertainties. For a 
lower assumed temperature for the foreground cluster, we fit 
a higher temperature for the background cluster. For example, 
if we assume that the temperature of the foreground cluster is 
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Fig. 2. The spectrum of the NE (left panel) and S W (right panel) X-ray peak extracted from circular region with a radius of 1 .5'. 
The continuous line represents the best fit model. The EPIC/pn and the coadded EPIC/MOS spectra are indicated. The position 
of the Fe K line in the spectrum from the NE peak shows that the emitting cluster is at a redshift of z = 0.44. 



Table 2. Fit results for the spectra extracted from the SW X- 
ray peak, from the apparent "tail" of the NE peak toward SW, 
and from the NE X-ray peak. The spectra extracted from the 
SW peak and from the "tail" region were fitted with a ther- 
mal model with th e A^h fixed to the Galact ic value (A^h 



1.47 x 10 20 cm' 2 , Dicke y & Lockmanlll990l) . The NE X-ray 



peak was fitted with two thermal models: one for the emission 
from A3 128 in the foreground and one for the background clus- 
ter. In this fit the redshift, temperature, and the Fe abundance of 
the foreground cluster were fixed to the best fit values from the 
"tail" region. The reported best fit parameter of the NE peak 
describe the properties of the background cluster. The temper- 
ature is given in keV, the emission measure Y is given in units 
of 10 65 cm -3 and the Fe abundance i s given with respect to the 
proto-solar values of lLoddersI (120031) . 



Par. 


SW 


"tail" 


NE 


^2=0.06 


9.72 + 0.10 


6.80 ± 0.07 


3.4 + 0.4 


Y z =0M 






303 ± 34 1 ' 


kT 


3.36 + 0.05 


3.50 ± 0.08 


5.14 + 0.15 1 


Fe 


0.69 ± 0.04 


0.49 ± 0.04 


0.47 ± 0.08+ 


z 


0.060 ± 0.001 


0.058 ± 0.003 


0.444 ± 0.003' 


X 2 /v 


652/501 


435/398 


943/800 



values determined for the background cluster. 



3.0 keV instead of 3.5 keV, the best fit temperature of the back- 
ground cluster will be 5.8 keV instead of 5.1 keV. 

We note, that there is no evidence of emission from the 
background cluster in the "tail" region. From the spectral anal- 
ysis we find that its contribution in this region is less than 1%. 



3.1.2. Radial profiles 

We determine the emission measure, temperature, and Fe abun- 
dance profiles for both X-ray peaks. We extract the spectra 
from circular annuli with outer radii of 0.5', 1.5', 2.5', and 3.5'. 



Table 3. Radial profiles determined for the S W X-ray peak. The 
temperature is given in keV, the emission measure Y is given in 
units of 10 65 cm 3 , and the Fe abundance i s given with respect 
to the proto-solar values of lLoddersI (12003b . 



Par. 


0.0' -0.5' 


0.5'-1.5' 


1.5'-2.5' 


2.5'-3.5' 


Y 

kT 

Fe 


1.99 + 0.05 
3.36 + 0.11 
0.87 + 0.10 


7.91+0.09 
3.36 + 0.06 
0.67 ± 0.04 


8.67+0.10 
3.17 + 0.06 
0.48 ± 0.06 


9.57 + 0.12 
3.38 + 0.07 
0.46 ± 0.04 


X 1 v 


269/229 


1125/959 


1183/1072 


1273/1137 



Table 4. Radial profiles determined for the NE X-ray peak. The 
metallicity and the temperature of the foreground cluster were 
fixed to 0.49 solar and 3.5 keV, respectively. The emission mea- 
sure Y is given in units of 10 65 ctrT 3 and the temperature is 
given in keV. The F e abundance is g iven with respect to the 
proto-solar values of lLoddersI (12003 ). The f := oM /./j=o.06 indi- 
cates the ratio of the fluxes from the background cluster and 
A3 128 in the given extraction region. The superscript "/" indi- 
cates that the parameter value was fixed during the fitting. 



Par. 


0.0' -0.5' 


0.5'-1.5' 


1.5'-2.5' 


2.5'-3.5' 


Y z =0M 


0.56 ± 0.20 


2.7 ± 0.4 


4.7 ± 0.4 


5.6 + 0.3 


Y Z =0A4 


72+ 15 


178 ± 27 


154 ± 34 


< 140 


kT 


5.0 + 0.3 


5.7 + 0.3 


3.5 + 0.3 


¥ 


Fe 


0.56 + 0.18 


0.49 ± 0.09 


0.21+0.11 


Q.2 f 


/:=0.44 
/—0.06 


2.1 


1.2 


0.55 


0.09 




260/266 


827/749 


986/841 


1109/870 



In Table [3] we show the best fit parameters determined for 
the SW X-ray peak. The temperature profile of this peak is 
flat and there is no evidence for a cool core. However, the Fe 
abundance has a strong peak in the centre of the X-ray emis- 
sion, which coincides with an elliptical galaxy surrounded by a 
group of galaxies. The abundance peaks within the central 1.5', 
which corresponds to ~ 120 kpc. 
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Fitting a beta model ( Cavaliere & Fusco-Femianolll978 ) to 
the radial surface brightness distribution of the SW X-ray peak 
we find P _ 0.3 and a core radius of r c = 30 kpc. As already 
noted by Rose et al.l d2002l) this core radius is about an order 
of magnitude lower than the more typical values of ~250 kpc 
that are observed for non-cooling core clusters of galaxies. The 
P value is also low compared to the typical value of ~0.7. The 
small value of f3 shows that the surface brightness distribution 
of the core is broader than that usually observed for clusters, 
which indicates that we might see the emission of a group su- 
perimposed on the diffuse cluster emission. 

Using the parameters of the /3 model fit, and a global tem- 
perature value of kT = 3.4 keV, we find that the mass enclosed 
within the radius of 120 kpc is 1 .2x 10 13 M . The total gas mass 
within the same volume, determined using a central electron 
density of 1.2 x 10~ 2 cirT 3 , is ~1.3 x 10 12 M Q . Using the same 
parameters, the estimated total mass within a radius of 1 Mpc is 
al.l x 10 14 M Q . A radius of 1 M pc is approximately th e virial 
radius of a 3.4 keV cluster (e.g. lFinoguenov et al 
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Fig. 3. Spectrum of the galaxy associated with the radio source 
SUMSS J033057-52281 1 in the centre of the NE X-ray peak. 
On the y-axis we plot arbitrary count units (ADU). Lines, used 
for redshift determination, are shown with arrows. 



The total Fe mass enclosed within the radius of 120 kpc, the 
region within which the Fe abundance peaks, is ~1 .5 x 10 9 M Q . 
Assuming a fiat Fe abundance distribution in A3 128 of 0.45 
solar, the excess Fe mass in the SW X-ray peak is ^5xl0 8 M Q . 



The radial profiles determined for the background cluster 
seen as the NE X-ray surface brightness peak are shown in 
Table [4] For the temperature and Fe abundance of the fore- 
ground cluster we assume 3.5 keV and 0.49 solar, respectively. 
The free parameters in the fit are the temperature and the Fe 
abundance of the background cluster, and the emission mea- 
sures of both clusters. We see an indication of a temperature 
drop in the inner 170 kpc, then the temperature peaks, and out- 
side 500 kpc it drops again. The spectra indicate that the Fe 
abundance peaks in the cluster core. In Table |4] we also show 
for each annulus the ratio of the flux of the background cluster 
with respect to the flux of the foreground cluster. In the inner- 
most extraction region the emission of the background cluster 
clearly dominates, between 0.5'-1.5' the contributions of both 
clusters are similar. Between 1 .5'— 2.5' the emission of the fore- 
ground cluster dominates, and between 2.5'-3.5' we can only 
determine an upper limit for the emission of the background 
cluster. The total 0.3-10.0 keV luminosity of the z - 0.44 clus- 
ter within the radius of 3.5' is 6.9 x 10 44 ergs s . 

To obtain a rough estimate of the total mass of the back- 
ground cluster, we fit a beta model to the radial profile of its 
best fit emission measure per arcminute. The best-fit parame- 
ter values are r c = 157 kpc and /? = 0.41. Assuming a global 
temperature of 5.14 keV (the best fit global temperature from 
Sect. 13. 1 . lb we obtain a total mass within a radius of 1 .5 Mpc 



(estim ated rsoo of a 5.1 keV cluster, e.g. see iFinoguenov et al 
2001) of 3.4 x 10 14 M . We note that due to the uncertainties 
associated with the subtraction of the foreground cluster, and 
because the core of the background cluster is far from being re- 
laxed and it is not single-temperature, this derived mass is only 
a very rough estimate. 



Fig. 4. Combined R-band image of the galaxy associated with 
the radio source SUMSS J033057-52281 1 in the centre of the 
NE X-ray peak. The total exposure time is 360 s. Note the arc 
to the Southwest of the galaxy. The scale of 1' is indicated on 
the image. 



3.2. Optical observation of the background cluster 

Optical spectroscopic observations of the gala xy associated 



with the radio source SUMSS J033057-52281 1 dMauch et al 



120031) were done with the 6.5 m Magellan I Baade telescope, 
using the Inamori Magellan Areal Camera and Spectrograph 
(IMACS). We used the short f/2 camera with a 27' x 27' FOV 
and 0.2" pixel scale. In spectroscopic mode this setup with 300 
1/mm grism and 1 .2" slit gives a spectral resolution of a: 6 A. 

The spectrum obtained is shown in Fig. [3] It is typical for an 
early type galaxy and does not show detectable emission lines. 
Using the standard set of absorption lines, marked in Fig. [3] 
with arrows, the redshift z = 0.43961 + 0.00014 was derived. 
This value is in a very good agreement with the X-ray deter- 
mined redshift thus providing immediate confirmation that the 
observed X-ray emission is due to the emission of a distant 
cluster. 
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Table 5. The best fit values determined for the X-ray bright low 
entropy region and for the region at the centre of the galaxy 
distribution. The emission measure Y s per square arcminute 
is given in units of 10 64 crrT 3 arcmin" 2 and the temperatures 
are in keV. The Fe abundance is g iven with respect to the 
proto-solar values of lLoddersI 020031) . The entropy s is given 
in keV cm 2 and the pressure P is given in 10" 12 dyne cm" 2 . 



Par. 


X-ray bright reg. 


centre of gal. dist. 


Ys 


7.54 ± 0.08 


4.03 ± 0.06 


kT 


3.30 ± 0.08 


3.79 ±0.13 


Fe 


0.43 ± 0.04 


0.53 ± 0.07 


s 


223 


335 


P 


9.5 


7.3 


x 2 /y 


563/569 


448/460 



During the slit alignment procedure we also obtained a 
few direct R-band images of this field. The combined image 
is shown in Fig. [4] The field was photometric ally calibrate d 
using observations of Landolt standard stars dLandoitHl992l) . 
We measured the magnitude of the galaxy to be m R = 18.78 
(in 4"aperture), which is approximately what i s expected for a 



brigh test cluster cD galaxy at this redshift (e.g. Vik hlinin et al 
119981) . The absolute magnitude of the galaxy in R was esti- 
mated to be -23 .35 (adopting K and evolution corrections from 
Poggianti 1997 ). For comparison the absolute magnitudes of 
M 87 and Cygnus A are approximately -23.2 and -23.3 respec- 
tively - close to the absolute magnitude of SUMSS J033057- 
522811. 

The gravitational lensing arc is clearly seen around 
this galaxy (Fig. |4) at a radial distance of ~ 6.2". The 
enclosed mass within this radius can be estimated (e.g. 
iNarayan & Bartermannf l996) as 



M = 1.1 X 10 14 M C 



{— ) 2 l— ^ 
1 30"/ UGi 



Gpc/ 



(1) 



where 9 = 6.17" is the arc radius and D - D^D^jD^ is the 
combination of (angular diameter) distances from the observer 
to the lens and to the source D s , and from the lens to the 
source D^. In our case only = 1.17 Gpc is known. For 
D — 1 Gpc we obtain an upper limit to the total mass of M < 
5 x 10 12 M within -35 kpc of SUMSS J033057-522811.This 
value is about a factor of 2 h igher than the total mas s (within 



2002), derived 



a similar distance) of M 87 (Ma tsushita et al. 
using X-ray data. 

One can therefore conclude that optical data strongly 
suggest that SUMSS J033057-522811 is a massive elliptical 
galaxy with the parameters characteristic for most bright cD 
galaxies in the local Universe. 

In the Magellan image we also see an excess of the surface 
density of faint galaxies near this cD galaxy, most of which are 
probably members of the distant cluster. 

3.3. 2D maps of thermodynamic properties 

For making 2D maps of spectral parameters, we select regions 
according to the cluster surface brightness in the 3.5-7.5 keV 



band. Considering a optically thin plasma model iForman et al 
(2007) showed that for Chandra ACIS-I spectral response the 
flux F in the 3.5-7.5 keV band coming from a unit volume 
with a given pressure depends only weakly on the gas temper- 
ature (over the 1-3 keV temperature range and for metallicity 
of -0.7 solar). For XMM-Newton EPIC and for the range of 
temperatures characteristic for A3 128 (3-5 keV) this is also 
approximately correct: a factor tj relating the pressure and the 
square root of the flux 77 = (e(T)/T 2 ) 1 ^ 2 decreases monotoni- 
cally as a function of temperature by about 20% (for metallic- 
ity of 0.45 solar). Here e(T) is the gas emissivity in the 3.5-7.5 
keV band (with account for the spectral response of the instru- 
ment). Therefore, areas with similar surface brightness in the 
3.5-7.5 keV band, provide us with contiguous regions with- 
out strong pressure discontinuities. We note here that our final 
maps of spectral parameters depend very weakly on the proce- 
dure of region selection. 

We then use the Voronoi tessellatio n method 
dCappellari & Copinl 120031 iDiehl & Staffer! 120061) to further 
bin the selected areas with a stable minimum signal-to-noise 
ratio of S /N = 33, which is needed for relatively accurate 
temperature determination. We identify 213 independent 
regions from which we extract the spectra. For each region 
we compute a spectral redistribution file and an ancillary 
response file. We fit the spectrum of each bin individually with 
a single temperature thermal plasma model. The abundances 
of all elements except Fe in our model are fixed to 0.4 times 
the solar value. The emission measure, temperature, and Fe 
abundance are free parameters in the fit. From the best fit 
emission measure and temperature, we calculate for each 
extraction region the density «, entropy s = kT/n 2 ^ 3 , and the 
pressure P = nkT, assuming a length scale along our line of 
sight of 1 Mpc. As a first step, we use a constant length scale in 
order to calculate the projected values of the thermodynamic 
properties without making any assumptions about the cluster 
centre and gas distribution. 

In order to subtract the emission of the background clus- 
ter from the 2D maps, we include in the fitted model for each 
bin with a mean distance smaller than 3.5' from the NE X-ray 
surface brightness peak, an additional thermal component. We 
fix its temperature, metallicity, and redshift to the mean val- 
ues determined for the background cluster from the global fit: 
5.14 keV and 0.47 solar, respectively (see Sect. 13.1. Ti l. Using 
the best fit beta model to the background cluster (see Sect. 
13.1.2b we calculate the normalization of the emission of the 
background cluster for each bin. By including this component 
in the fitted model, we effectively subtract its emission and 
determine the best fit parameters for the foreground cluster. 
However, since the core of the background cluster is not ra- 
dially symmetric and it has a strongly elongated morphology, 
at radii smaller than 0.8' from the NE peak its emission cannot 
be properly subtracted. Therefore, we exclude the bins with a 
distance smaller than 0.8' from the NE peak from the analysis. 

At the top left panel of Fig. |5]we show the map of the pro- 
jected emission measure per square arcminute integrated along 
the line of sight, after the emission of the background cluster 
was subtracted. We see, that the surface brightness of the X-ray 
emission of A3 128 has a peak at about 1.5' West-Southwest 
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Fig. 5. Top left: The projected emission measure per arcminute, in units of 10 64 cirT 3 arcmin~ 2 , integrated along the line of sight, 
after the emission of the background cluster was subtracted. Top right: Temperature map derived from spectral fitting in units of 
keV. Bottom left: Map of the pressure derived from the best fit temperature and emission measure, assuming a length scale along 
our line of sight of 1 Mpc. The units are 1CT 12 dyne cirT 2 . Bottom right: Map of the entropy derived from the best fit temperature 
and emission measure, assuming a length scale along our line of sight of 1 Mpc. The units are keV cm 2 . The white and black 
circles indicate the high surface brightness low entropy region and the region at the centre of the galaxy distribution, respectively. 
These we use for subsequent spectral analysis. The X-ray isophotes from the wavelet decomposed image are overplotted. 
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Fig. 6. Left panel: Density map of the ICM. Central panel: Map of the entropy of the ICM. Right panel: The pressure map of 
the ICM. The ther mod ynamic propert i es are calculated using volume estimates for the spectral extraction regions described in 
Henr y et al (2004) and Mahda vi et alJJ2005h . The volumes are calculated assuming two centres, one at the SW X-ray peak, and 
one at the X-ray bright low entropy region at the possible centre of the gravitational potential of A3 128. 
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from the original NE peak. This region, as we show in the 
top right panel of Fig. [5] has a temperature between 3.0-3.5 
keV. While the pressure is relatively high in this region (see 
the lower left panel of Fig. [5j, the entropy is low (lower right 
panel of Fig. [5]). This indicates that this region may be the cen- 
tre of the gravitational potential of A3 128. The centre of the 
galaxy distribution of A3 128 is abou t 2.8' South of this region 
dPalton et al.lll997tlRose et~ai1l2002|) . 



In order to better compare the spectral properties of these 
two regions: the X-ray bright low entropy region and the region 
at the centre of the galaxy distribution, we extract a spectrum 
from a circular area with a radius of 1 .26' centred on the new 
X-ray brightness peak (or = 3 h 30 m 40 s , 5 = -52°28'50") and a 
spectrum from a circular area with a radius of 1.5' centred on 
th e optical centre (a = 3 h 30 m 43 s , 6 = -52°31'30") reported 
bv lRose et al.1 020021) . The spectral extraction regions are indi- 
cated by circles on the maps in Fig. [5] The best-fit values of 
the single temperature thermal fits to the spectra are shown in 
Table [5] We also show in the table the entropy and pressure 
values calculated for the volume of a sphere with a radius cor- 
responding to the radius of the extraction region. The results 
confirm that the high surface brightness region has a factor of 
~1.5 lower entropy than the area at the optical centre of the 
galaxy distribution. The pressure at the high surface brightness 
region is a factor of ~1.3 higher than the pressure at the cen- 
tre of the galaxy distribution. The fit results indicate that the 
Fe abundance of the X-ray bright region is lower than the Fe 
abundance at the centre of the galaxy distribution. Under the 
assumption of hydrostatic equilibrium, the pressure peak com- 
bined with low entropy are unambiguous signs of the centre of 
the gravitational potential. 



The X-ray bright region at the possible centre of the gravita- 
tional potential of A3 128 in the NE is separated from the bright 
SW core by a surface brightness depression. The pressure map 
indicates that the dark matter potential well associated with the 
S W core is connected with the dark matter potential in the NE. 
This is possibly the region where the dark matter potentials of 
the SW core and of A3 128 overlap. We note that the "bridge" 
between the SW and NE in the pressure map coincides with an 
apparent chain of galaxies seen in the Digitalised Sky Survey 
image. Unfortunately, non of these galaxies has a known radial 
velocity. 



In Fig.|6l we show density, entropy, and pressure maps de- 
termined assuming the X-ray emitting gas is associated with 
two gravitational potentials. The thermodynamic properties are 
calculated using volum e estimates for th e spe ctral extraction 



region s described in iHenry et alJ (120041) and iMahdavi et al 



(2005). The volumes are calculated assuming two centres, one 
at the SW X-ray peak, and one at the X-ray bright low entropy 
region at the possible centre of the gravitational potential of 
A3 128. For each polygon the closest of the two centres is as- 
sumed to be the centre in the volume calculation. 



4. Discussion 

4.1. The background cluster associated with the NE 
X-ray peak 

While at the expected energy we do not detect the strong Fe K 
line in the spectrum extracted from the NE X-ray peak, we de- 
tect line emission at the energy corresponding to the line energy 
of Fe K emission redshifted by z = 0.44. Therefore we con- 
clude that the NE X-ray peak observed toward A3 128 is not as- 
sociated with the surviving ICM of a group falling in super soni- 
cally to the cluster as previously thought ( Rose et al]|2002 ), but 
it is a distant luminous cluster of galaxies at redshift z = 0.44. 
Subsequent optical spectroscopic observation of the distant ra- 
dio bright galaxy in the centre of the NE X-ray peak (radio 
source SUMSS J033057-52281 1) with the Magellan telescope 
revealed a redshift of z = 0.44 confirming its association with 
the cluster seen in X-rays. 

The red magnitude of the galaxy otr = 18.5 is typica l 



for cD galaxies at this redshift (e.g. Vikhli nin et alJ [1998). 
The Magellan observation also confirmed the gravitational arc 
around the galaxy, the presence of which was previously sug- 
gested b y the images obtained with the 0.9 m CTIO tele- 
scope dRose et al.l l2002). We estimate the enclosed total mass 
within the radius of the lens to be < 5 x 10 12 M , which is 
within a factor of 2 consistent with the total mass of M 87 
at a similar radius. The optical data thus strongly suggest that 
SUMSS J033057-52281 1 with properties characteristic for cD 
galaxies is the dominant galaxy of the cluster in the back- 
ground. 

The X-ray morphology of the background cluster is 
strongly elongated showing that its core is not relaxed and 
that it is possibly undergoing a merger. However, its observed 
X-ray luminosity of L = 6.9 x 10 44 erg s _1 , temperature 
kT = 5.14 keV, and derived mass M = 3.4 x 10 14 M agree 
well with th e generally observed Lx-T and M—T scaling re 



lations (e.g. IWu et al. I ll999t IPratll 12006). Further deep optical 
observations are needed to characterise the properties of the 
background cluster. 

4.2. The nature of the SW surface brightness peak 

The S W surface brightness peak is centred on an apparent com- 
pact group of galaxies. Compared to other clusters of galaxies 
with similar temperatures, its core radius is small. Such small 
core radii are observed only in clusters with cooling cores, or 
in clusters undergoing merging events where the cluster pro- 
file is truncated. While the radial Fe abundance distribution 
has a strong peak at the centre of the group, the temperature 
distribution is flat, with no indication for a cool core. The 
cooling time in the centre of the SW core, assuming the gas 
cools isobarically, is 4.5 x 10 9 yr. Its estimated total mass of 
M500 = 1.1x10 14 is a fact or of 2 lowe r than the expected mass 
of a 3.4 keV cluster (e.g. iPratt 120061) . But since the value of 
/} = 0.3 was determined by fitting the surface brightness pro- 
file at radii much smaller than rsoo, the estimated M500 is highly 
uncertain. The radial velocity of the dominant gal axy of this ap- 
parent group (ENACS 75) is V = 19252 km s~ 2 dKatgert et al.1 
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19981) . which is higher by » 1500 km s than the mean radial 
velocity of galaxies in A3 128. Unfortunately, we know the ra- 
dial velocity for only one more galaxy of the apparent group 
(ENACS 78), which is lower than that of t he dominant galaxy: 
= 18380 km s - 



V 



(Katg ertet all 1 1998), but still higher by 



ss 500 km s _1 than the mean radial velocity of the cluster. The 
sound speed corresponding to the mean cluster temperature of 
~3.5 keV is 960 km s _1 , which means that if the SW peak is 
the remaining core of a group or a cluster merging with A3 128 
at a relative velocity of « 1500 km s _1 , then the merger is su- 
personic. 

On the smaller spatial scales, the X-ray emission of the 
SW peak appears to follow the galaxy mass distribution as 
it is slightly displaced toward the West and clearly peaks on 
the galaxy ENACS 75. However, the association of the galaxy 
ENACS 75 with the SW X-ray peak in radial velocity is diffi- 
cult. The redshift of the SW peak determined from the energy 
centroid of the Fe K line corresponds to a radial velocity of 
V = 18000 + 300 km s -1 , and the redshift of the "tail" region, 
which we might consider as the diffuse emission of A3 128 cor- 
responds to a radial velocity of V = 17400 + 900 km s _1 . 
However, the systematic uncertainties on these values are larger 
than the quoted statistical errors. The EPIC detectors are known 
to have gain problems which are both time and position depen- 
dent. The absolute error in the redshift determination can be as 
large as 1500 km s (Simionescu et al. in prep.). Therefore, 
the radial velocity of the hot gas is within the systematic uncer- 
tainties consistent with the radial velocity of the galaxies. 

If the SW peak is a core of a group or of a cluster that 
moves through the ICM of A3 128, than its gas is being stripped 
by ram pressure. The gas of the infalling cluster is stripped at 
radii where the thermal pressure of the hot gas in the infalling 
group is too small to balance the sum of the thermal and ram 
pressure of the cluster ICM. This happens at radii larger than 
the radius where the thermal pressure of the infalling cluster 
is equal to the thermal press ure at the stagnation point (e.g. 
Markevitch & Vikhlininll2007h . The ratio of thermal pressures 
at the stagnation point po and in the free stream, p\, for Mach 
numbers relative to the sound spe ed in the free stream region 
M > 1 is (Land au & Lifshitzfr9 59): 



Pi 



fy + 1 



M \y — 



y-i 

2M 2 



(2) 



where y = 5/3 is the adiabatic index of the gas. For the pressure 
in the free stream region, pi, we assume the value determined 
assuming a density of n — 1 x 10~ 3 cm -3 and a temperature of 
3.5 keV. For the pressure of the SW core at the given radius, 
we assume a density profile of n{r) = hq{\ + (r/r c ) 2 y^, where 
n = 1.2 X 10~ 2 cnr 3 , r c = 30 kpc, and f3 = 0.3. We find that at 
a Mach number of M — 1.5 (inferred from the radial velocity 
difference of ~ 1500 km s" 1 between ENACS 75 and A3 128, 
and from the sound speed v s = 960 km s in the ICM), the 
gas outside the radius of r = 90 kpc should be stripped, and 
only the ICM of the merging core inside of this radius should 
be surviving. 

Such a merger would produce a tail of low entropy and 
high Fe abundance. However, in the images and in the maps 



of the thermodynamic properties we do not observe any ob- 
vious signs of a merger for the SW peak. Moreover, we ver- 
ified that the Fe abundance distribution is consistent with be- 
ing symmetric and is not enhanced in any direction. These ob- 
served properties might indicate that the merger is occurring 
just along our line of sight. The merger is compressing and 
heating the ICM, which might explain the observed lack of a 
cool core. A shock with a Mach number M = 1.5 heats the 
gas by a factor of 1.5, which means that for the cluster tem- 
perature of 3.5 keV the gas temperature at the stagnation point 
should be 5.25 keV. If we fit the spectrum extracted from the 
circular region with a radius of 1.5' centred on the SW core 
with two thermal models, with the temperature of one thermal 
component fixed to kT = 5.25 keV, the fit improves compared 
to the single temperature model (see Sect. 13.1. Ti l. The reduced 
X 2 improves from 1.30 to 1.16. We obtain a best fit tempera- 
ture of kT = 2.23 + 0.17 keV for the cooler component, and 
best emission measures of Y = (5.1 + 0.5) x 10 65 cirT 3 and 
Y = (4.7 + 0.5) x 10 65 cirT 3 for the hotter and cooler gas, re- 
spectively. This indicates that much of the emission of the SW 
core might be due to the shocked gas from a merger happening 
along our line of sight. As the enriched stripped gas is trailing 
behind the core we are seeing it in projection which explains 
the enhanced Fe abundance within the projected distance of 
-120 kpc. The velocity difference between the two galaxies of 
the compact group (ENACS 75, ENACS 78) may also be ex- 
plained by tidal stretching of the group along our line of sight 
during the passage of the group through the cluster. 

Alternatively, the SW core might be a group still falling 
toward A3 128, which we observe projected on the ICM of the 
cluster. In this case the ICM of the group still did not start to 
interact with the ICM of the cluster, and ENACS 78 and 75 
are only projected close to each other by chance. However, this 
scenario does not explain the lack of the temperature gradient 
and of a cool core in the group. 



4.3. The diffuse ICM 

After subtracting the emission of the z = 0.44 cluster, we iden- 
tified a new region with an enhanced surface brightness to the 
West-Southwest of the centre of the background cluster. This 
region has also an enhanced pressure and low entropy. It is at 
a distance of ~2.8' from the centre of the galaxy distribution 



given by lRose et al.l (120021) . However, the entropy of the ICM at 
the centre of the galaxy distribution is a factor of 1 .4 higher and 
the pressure is lower than that of the region with the enhanced 
X-ray surface brightness. Because of the great dynamical com- 
plexity of the system, with infalling groups and filaments iden- 
tified in position-position and position-redshift diagrams, the 
determination of the cluster centre from the galaxy distribution 
is highly uncertain. Based on the thermodynamic properties of 
the ICM, we conclude that the enhanced surface brighness re- 
gion, centred at a = 3 h 30 m 40 s , 5 = -52°28'50", is the centre 
of the gravitational potential of the cluster A3 128. The posi- 
tion of this region is also more consistent with the centre of the 
extended low surface brightness X-ray emission. 
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The images and the maps of thermodynamic properties do 
not reveal any obvious shocks. We only see two candidates, 
one to the Northeast of the background cluster and one to the 
Southwest of SW core. Both shock candidates exhibit temper- 
ature drops by ~1 keV associated with a drop in entropy, and 
a factor of 2 drop in pressure. However, the surface brightness 
in these two areas is low and we cannot confirm the presence 
of surface brightness discontinuities in the regions where the 
maps indicate the p resence of the shocks. 

As proposed bv lCaldwell & Rose! (11997b . A3 128 may have 
had encou ntered a merger w ith A3 125 in the past, and as dis- 
cussed by iRose et al. (2002) the optical redshift data reveal a 
number of groups, some of them tidally distended into fila- 
ments after a close passage through A3 128. The X-ray data 
combined with the optical redshifts suggest an ongoing merger 
with a group. Moreover, the inconsistency between the clus- 
ter centre determined based on the distribution of galaxies and 
based on the thermodynamic properties of the ICM, and the 
lack of bright galaxies at the newly identified centre of grav- 
itational potential, further highlight the view that A3 128 is a 
dynamically young, unrelaxed system. The unrelaxed nature of 
A3 128 can be attributed to its location in the high density envi- 
ronment of the Horologium-Reticulum supercluster. 



5. Conclusions 

We have analyzed new deep XMM-Newton EPIC data of the 
cluster of galaxies Abell 3128 located in the Horologium- 
Reticulum supercluster. We found that: 

- The Northeast X-ray peak observed toward A3 128 is a dis- 
tant luminous cluster of galaxies at redshift z = 0.44. 

- The properties of the distant radio bright galaxy in the cen- 
tre of the NE X-ray peak indicate that it is the cD galaxy of 
the cluster in the background. We detect a gravitational arc 
around the galaxy. 

- The properties of the Southwest X-ray peak suggest that it 
is the core of a group merging with A3 128 along our line 
of sight. 

- Based on 2D maps of thermodynamic properties of the 
ICM determined after subtracting a model for the back- 
ground cluster, we conclude that the enhanced surface 
brightness region at a distance of ~2.8' from the centre 
of the galaxy distribution is the centre of the gravitational 
potential of the cluster A3 128. The inconsistency between 
the cluster centres determined based on the distribution of 
galaxies and based on the thermodynamic properties of the 
ICM further highlights the view that A3 128 is a dynami- 
cally young, unrelaxed system. 
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